Interaction effects in the transport of two-dimensional holes in GaAs 
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The power-law increase of the conductivity with temperature in the nominally insulating regime, 
recently reported for the dilute two-dimensional holes cond-mat/0603053 , is found to systematically 
vary with the carrier density. Based on the results from four different GaAs heterojunction- insulated- 
gate field-effect-transistor samples, it is shown that the power law exponent depends on a single 
dimensionless parameter, the ratio between the mean carrier separation and the distance to the 
metallic gate that screens the Coulomb interaction. This dependence suggests that the carriers form 
a correlated state in which the interaction effects play a significant role in the transport properties. 

PACS numbers: 73.40.-c, 73.20.Qt, 71.27.+a 



Charge transport in two-dimensional (2D) electron sys- 
tems [l| provides a unique means of studying the inter- 
play between disorder and electron-electron interactions. 
This problem is fundamentally important and remains a 
subject of intense investigation Q. While noninteracting 
2D electrons are generally believed to form the Anderson 
insulator 0, 0] , the situation appears to be much more 
complex in the presence of interactions 0, IE @ • -^ e " 
cent theoretical studies 0, ^| emphasize the importance 
of collective phenomena around the so-called metal-to- 
insulator transition (MIT). 

In practice, since the ratio of the Coulomb interaction 
to the kinetic energy increases with lowering the den- 
sity, samples with most dilute carriers are best suited 
for probing the collective phenomena. However, reduc- 
ing the carrier density n runs into a risk of increasing 
the carrier separation oc n -1 / 2 beyond the single-particle 
localization length £, in which case the interaction effects 
become overshadowed by the single-particle localization. 
Thus, a sufficiently clean 2D environment is another re- 
quirement that has to be met to uncover the underlying 
interaction effects. 

The experimental progress on studying the transport 
of 2D systems has been greatly influenced by the sam- 
ple quality. Very early experimental results in Si-devices 
demonstrated the activated transport consistent with 
the Anderson insulator Q: the Arrhenius conductivity 
a ~ e~ TA l T at high temperatures, and the softer tem- 
perature dependence, a ~ e~' T / T ' at lower tempera- 
tures, with v = 1/3 corresponding to the variable-range 
hopping (VRH) scen ario [ll|], and v = 1/2 to the effect 
of the Coulomb gap 0. However, in the mid-1990's, 
experiments performed in much cleaner 2D electrons in 
Si-MOSFETs showed both metal-like (da/dT < 0) and 
insulator-like (da/dT > 0) conductivity behavior, de- 
pending on whether the density is above or below a cer- 
tain critical value n c ^3] . Although the existence of the 
metallic regime at T — » is still debated, the transport 
on the insulating side generally remains activated, in ac- 



cord with the Anderson localization picture. 

The experimental findings on the insulating side of 
the MIT has radically changed with the adoption of 
the undoped GaAs/AlGaAs heterojunction-insulated- 
gate field-effect transistors (HIGFETs). Recent exper- 
iments |l4l 1 151 demonstrated that the conductivity can 
become non-activated, while preserving the "insulating" 
sign da/dT > 0. In 2003, a close-to-linear dependence 
a oc T was first observed in 2D holes in a p-GaAs 
HIGFET device for densities down to p = 1.5 x 10 9 cm~ 2 
|l4| . Subsequent experiments in the devices of the same 
kind not only confirmed this observation for similar car- 
rier densities, but also revealed a more general power- 
law- like temperature dependence, a oc T a , with a vary- 
ing e xpon ent 1 % a < 2, at sufficiently low tempera- 
tures [l5j. Remarkably, such a behavior persists even 
for a record- low density of p — 7 x 10 8 cm -2 , in which 
case the Coulomb energy is about 100 times greater than 
the nominal Fermi energy, and the Fermi wavelength 
Xf = \/2n /p — 0-95 fj,m approaches a macroscopic scale. 

In this Letter we present a comprehensive study of this 
surprising power law T-dependence based on data col- 
lected from four different p-typc HIGFET samples that 
only differ by the structural barrier thickness (the dis- 
tance d from the 2D layer to the metal gate, Table [i}. 
The procedures of sample preparation and the measure- 
ment details are described in Ref. | 15| for the first three 
samples, while the data from the fourth sample is drawn 
from Ref. 0] for comparison. The fitting of the conduc- 
tivity for the temperatures 35 < T < 200 mK to 

a/a Q = Go + (T/T Q ) a , a Q = e 2 /(2irh) (1) 

yields a sample-dependent exponent a(p) that grows with 
decreasing hole density p, while the T-independent term 
remains negligible, Go<7q <C a(T). Moreover, we find 
that a systemically depends on the single dimensionless 
parameter k = a/d, where a = (-np)^ 1 / 2 is the Wigner- 
Seitz radius (the mean carrier distance is about 2a). 
Such a single-parameter dependence, tied to the screen- 
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FIG. 1: Conductivity temperature dependence in the log-log 
scale for a set of hole densities, for the samples (a) #2, (b) 
#3, (c) #4, and (d) the sample from Ref. |l4j|. The scattered 
points in (b) are the dc results for p — 8 x 10 8 cm -2 . 



ing length ~ 2d for the Coulomb interaction rather than 
to the sample-specific localization length, strongly indi- 
cates an important role played by the electron interac- 
tions. 

The aim of this work is to analyze the power-law-like 
temperature dependence <j(T) and to examine the rele- 
vance of the interaction effects. In Fig. ^ we show the 
log-log plots of <r(T) obtained from four different sam- 
ples: (a) #2 and (b) #3 are samples from the same wafer, 
cooled down to the lowest temperature of around 80mK; 
(c) sample #4 is cooled down to 35 mK; and (d) is the 
data from Ref. 0] with the lowest T of 65 mK. Each 
panel in Fig. ^ contains curves for a number of hole den- 
sities. For the density above critical, p c ~4x 10 9 cm~ 2 , 
the transport is metal-like Below we focus on the 



TABLE I: The barrier thickness d (distance to the metallic 
gate) for four different samples 



Samples 


#2 


#3 


#4 


Noh et al. 


Barrier d 


600 nm 


600 nm 


250 nm 


500 nm 



opposite, low-density "insulating" regime. 

There are three features common to all four sam- 
ples, (i) Although the sign da/dT > 0, the conductiv- 
ity never becomes activated |l5J because an activated In- 
dependence in a log-log plot would have shown a strong 
downward bending for T below the activation tempera- 
ture, (ii) For low densities p < 2 x 10 9 cm -2 , the slopes 
d log a/d log T are roughly temperature- independent at 
temperatures ^ 200 mK and increase with decreasing 
density, (iii) Finally, although the non-activated con- 
ductivity is a signature of the extended states, the con- 
ductivity values are 1-2 orders below e 2 /(2Trh) (the 2D 
analog of the minimal metallic conductivity 16] ). 

In what follows we investigate the density dependence 
of the slope a — d log a/d log T by fitting a{T) to the 
simple formula JQ) for the low density curves, p < 3 x 
10 9 cm -2 . The fitting parameters Go and Tq are plotted 
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FIG. 2: Density dependence of the fitting parameters Go and 
To from Eq. Ql 
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in Fig. [5] The T-independent term Go remains around 
zero up to p ~ 2 x 10 9 cm -2 for the samples #2, #3, and 
#4, indicating a fairly good power law dependence: a oc 
T a . On the other hand, for the sample from Ref. 0], the 
constant Go is more significant. We stress that, although 
the high-temperature behavior in the sample from Noh 
et al. looks quite linear, a ~ A + BT, as reported in 
Ref. 0| > at T < 200 mK the conductivity curve exhibits 
a noticeable bending at T < 200 mK as demonstrated in 
the linear scale plot in the inset of Fig^(d), and is best 
fitted to Eq. (Q) with a ^ 1 [see Fig. |3] . The parameter 
To for all four samples has a trend of a slow decrease with 
increasing density from about 400 mK to about 300 mK 
for densities up to p ~ 2 x 10 9 cm~ 2 . 

Fig. |3 shows how the exponent a depends on the pa- 
rameter 2a, which is approximately the mean carrier 
spacing. The results from samples #2 and #3 fall ap- 
proximately onto a single curve in which a varies from 
1.1 to about 1.5. The results from both samples #4 and 
Ref. qualitatively follow the same trend albeit the 
values of a are greater by about 0.6 and 0.1 respectively. 
This shift in a(p) relative to that for samples #2 and 
#3 motivates us to look into the structural differences 
between the samples, which is primarily in the barrier 
thickness d. As shown in Table [i] samples #2 and #3 
have the same d, while the values of d for the other two 
samples are notably different. Furthermore, the values of 
a increase with the decrease in d for a given density. This 
trend points at the role of the screening by the metallic 
gate. 

In HIGFETs, the metallic gate at the distance d from 
the 2D hole layer screens the 1/r interaction down to 
1/r 3 when r > 2d. For the lowest densities, the interac- 
tion becomes effectively short-ranged, with the relevant 
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FIG. 3: Density dependence of the exponent a for all the 
samples. The dotted lines are guides for the eye. 
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FIG. 4: The exponent a as a function of the ratio k = a/d. 
Inset: a 1 ^ as a function of T/To(p). 



parameter being the ratio k = a/d between the carrier 
spacing w 2a and the screening radius 2d. For our mea- 
surements, this ratio can be continuously varied in the 
range 0.1 < k < 0.85; the sample of Noh et al. also falls 
into this range. The effect of the gate screening becomes 
apparent in Fig. 01 where all the power law exponents a 
are plotted as a function of k. Remarkably a(n) from all 
four samples fall onto a single curve within reasonable 
error bars. This curve tends to saturate for a > 2, while 
it most strongly varies when a « 1. Thus the linear de- 
pendence a oc T reported in Ref. [l4T | is most probably 
a crossover into an entirely different transport regime, 
rather than a universal signature at low T. [In the inset 
we plot cr 1 /" as a function of T/To(p) to illustrate the 
relative insignificance of the constant Go in Eq. Q).] 

The strong sensitivity of the transport to the shape 
of the interaction potential is not entirely surprising. 
Indeed, the carriers at these densities are delocalized 
and are very strongly interacting (the r s value for p — 
1 x 10 9 cm -2 is 100 if one assumes the hole band mass 
m = 0.4m e ). A plausible way to think about such a 
system is by imagining a liquid or a strongly-interacting 
plasma. As it has been recently shown, the kinetic and 
thermodynamic properties of classical plasmas strongly 
depend on the same screening parameter k , although 
the classical arguments alone cannot explain the peculiar 
dependence Q . It is interesting whether the quantum ef- 
fects could manifest themselves at lower (so far inacces- 
sible) temperatures, in which case the system could be- 
come collectively localized Q , or whether the MIT takes 
place JIIj|- The dependence on n = a/d, where a and d 
can be independently controlled, suggests that, by vary- 
ing the density, one can continuously modify the state 
of the system. For sufficiently low densities (a > d), a 
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possibility of a reentry into the Fermi-Liquid (FL) was 
suggested 0. For larger densities, a sequence of mixed 
phases 0| was conjectured. 

We finally note that varying d not only changes the 
interaction range, but also the correlation length £dis of 
the disorder potential. The scattering off the surface im- 
perfections at the gate level is estimated to be negligibly 
small. In the case when the disorder is dominated by the 
charge impurities in the bulk, the gate screens the dis- 
order potential harmonics for length scales > d, so that 
£dis ~ d. In such a situation it becomes more difficult 
to differentiate the electron-electron interaction effects 
from those due to the change in the distribution of disor- 
der |2fJ,|2l|. The dependence of the conductivity on a/d, 
rather than on a and d separately, suggests that the elec- 
tron interactions are probably more important than the 
change in disorder. One natural possibility for the elec- 
tron interaction to enter is through the screening of the 
impurity field (such as via the RPA screening at smaller 
r s yl l22| ) j in which case the resulting effective disorder 
would depend on a/d as long as a < d. Another, equally 
plausible scenario is provided by assuming the hydrody- 
namic (viscous) flow of an electron liquid past the impu- 
rities whose size £dis > In this case the resistivit y is 
proportional to the viscosity of the 2D liquid [lH E^.l24j 
(that depends on k), and to the number of impurities, 
while the dependence on their shape and size (that may 
become affected by varying d) enters only under the log- 
arithm, according to the well-known Stokes paradox of a 
2D flow [13. 

In summary, by performing transport measurements 
on high quality 2D holes with densities down to 7 x 
10 8 cm -2 , we established the dependence of the power- 
law exponent a [Eq. Q)] on the ratio a/d between the 
Wigner-Seitz radius a and the distance to the metal gate 
d. We ascribe this dependence to the screening of the 
Coulomb interaction by the gate. We believe that our 
results provide direct evidence of the role of the electron- 
electron interaction in the 2D transport, and suggest that 
the transport is sensitive to the shape of the interaction 
potential controlled by the screening distance d. By vary- 
ing the ratio a/d, one can realize a strongly-correlated 
state of the 2D carriers with tunable properties. 
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